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Abstract 

In this paper we analyze a quartic tensor model with one interaction for a tensor of arbitrary rank. This model 
has a critical point where a continuous limit of infinitely refined random geometries is reached. We show that 
the critical point corresponds to a phase transition in the tensor model associated to a breaking of the unitary 
symmetry. We analyze the model in the two phases and prove that, in a double scaling limit, the symmetric 
phase corresponds to a theory of infinitely refined random surfaces, while the broken phase corresponds to a 
theory of infinitely refined random nodal surfaces. At leading order in the double scaling limit planar surfaces 
dominate in the symmetric phase, and planar nodal surfaces dominate in the broken phase. 


1 Introduction 

Tensor models mmmm generalize matrix models [5] and generate dynamical triangulations 00 in any dimension. 
The number of simplices in a triangulation diverges when the coupling constants are tunned to some critical values 
MESUn]. Sending at the same time the edge length of the individual simplices to zero, so as to keep the physical 
dimensionfull volume of the triangulation fixed, matrix and tensor models describe a theory of random infinitely 
refined spaces 00- 

Like matrix models E23, tensor models are endowed with a small parameter, 1/N (where N is the size of the 
tensor) and exhibit a 1 /N expansion [13, T4] [15] QTj, 4]- However, while in matrix models each order in 1/N 
corresponds to a fixed topology (i.e. to fixed genus surfaces) and the 1/N expansion is topological, the same can 
not hold true in higher dimensions. Indeed, for D > 3, there does not exist a single topological invariant that 
distinguishes topologies which could be used to index a topological expansion. Consequently the 1/N expansion is 
not topological in dimension higher than two. However, it is as close as one can get to one: the different topologies 
are explored progressively (i.e. at any order in 1/N only a finite number of topologies contribute) and at leading 
order only the spherical topology contributes. 

The triangulated spaces generated by a tensor model in rank D have topological dimension D. However, the 
emergent continuous geometry in the critical regime has typically a different effective dimension. While in the case 
of matrices the emergent geometry at leading order is a Brownian sphere mi Emi], with Hausdorff dimension 4 
and spectral dimension 2, for tensor models the emergent geometry at leading order is typically a branched polymer 
m (a continuous random tree in the mathematical literature), with Hausdorff dimension 2 and spectral dimension 
4/3. In order to uncover a genuine higher dimensional random geometry using tensor models one needs to go 
beyond the simple critical behavior at leading order. A recent proposition consists in boosting some subleading 
contributions [20] in tensor models. With these non canonical scalings one can write tensor models with a richer 
phase structure having notably a Brownian sphere phase. 

We follow here an alternative route towards a non branched polymer phase in higher dimension. Our approach 
is inspired by the double scaling limit mi mm mm mi- In the critical regime, the contributions of the lower 
orders in 1/N are enhanced. Taking at the same time the coupling constants to criticality and N to infinity one 
obtains a continuous phase in which many orders in 1/N contribute. In this way one can hope that the contributions 
of this larger family of graphs pile up and eventually lead to a new geometric phase. 

The double scaling limit of tensor models (D > 3) has so far been studied for the case of actions which are 
invariant under a global permutation of the indices of the tensor (i.e. for which all the indices of the tensor are 
treated on an equal footing). In this case the leading family in double scaling turns out to be exponentially bounded 
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for 3 < D < 6 and not exponentially bounded for D > 6 [24] (for matrices D = 2, the leading family in double scaling 
is not exponentially bounded). While not proven, it is very likely that all the leading order graphs in double scaling 
for 3 < D < 6 have the topology of a D dimensional sphere. However, as this family proliferates like a family of trees 
(the cherry trees [24]), it is very likely that for models possessing a color permutation symmetry the double scaling 
limit explored so far leads again to a branched polymer phase. 

In this paper we study in detail a tensor model with only one quartic interaction. The main difference with 
respect to the models studies previously is that in our case the indices of the tensor are not treated on an equal 
footing: as we will see below the action is not invariant under a global permutation of the indices of the tensor. 

We study the critical behavior of this model and show that (similarly to the color symmetric quartic model M) 
the critical point of the dynamical triangulation continuum limit corresponds to a phase transition in the tensor 
model. However, the similarities between the model with only one interaction and the color symmetric model end 
here. As the non symmetric model is simpler than the symmetric one, we are able to solve the equations of motion 
in both phases and identify the physical vacua of the theory on the two sides of the phase transition. We are then 
able to analyze the approach to criticality (in double scaling) on both sides of this phase transition. 

The main result of this paper is that, for the model with only one interaction, one obtains a genuinely new 
geometric phase in double scaling: instead of obtaining (at leading order in double scaling) a family of trees, one 
obtains (at leading order in double scaling) a family of planar graphs. The phase transition (in double scaling) is a 
phase transition between a phase of planar surfaces and a phase of planar nodal surfaces. 

Of course the emergent geometry obtained in double scaling is not yet genuinely D dimensional, but (very likely) 
is a Brownian sphere. However, this is a first step in the quest for a genuinely new random geometry in higher 
dimensions: we succeeded for the first time in finding a non branched polymer state of higher dimensional random 
geometry in a scaling regime of tensor models. 

Besides their importance for the invariant tensor models, our results could be relevant for group and tensor field 
theories Hi [Ml 1301133131133131 [Ml (SSJ. A regime parallel to the double scaling limit we discuss in this paper 
could be reached naturally by such models under the renormalization group flow [37] [38]. 

This paper is organized as follows. We introduce the quartic tensor model with one interaction in section El We 
obtain the vacua of this model and show that the critical point corresponds to a phase transition in section [3] We 
discuss the effective theory around a vacuum state in sections [4] and 0 Finally we analyze the random geometry 
underlying the effective theory in the two phases of the model in section [G] 


2 The model 


An introduction to tensor models in general and to the quartic models in particular can be found in mmmm- 
Let us consider a rank D tensor T a i a r>, where a c = 1.. .TV (and let use denote T a i a D its complex conjugated) 
transforming under the external tensor product of D fundamental representations of the unitary group U(N): 
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The position c e T> = {1,... D} of an index is called its color. A tensor model is a probability distribution for T 
and T invariant under the transformations in Eq. 0- The partition function of the quartic melonic tensor model 
symmetric under permutations of the colors is: 
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where, for any set C c 2?, we denote n c = (n c ,c e C ) and 5 n c n c = FlceC $n c n c - Each invariant V^ 4 ^(T, T) represents 
a gluing of four U-simplices, two positively and two negatively oriented. The coupling g counts the number of 
positively oriented simplices. When evaluating this partition function in perturbation theory one obtains Feynman 
graphs having (D + 1) colors which are dual to .D-dimensional triangulations [4] . 

This model can be studied in the intermediate field representation HE®]. Each quartic interaction can be 
represented as Gaussian integral over an auxiliary TV* x TV hermitian matrix: 
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normalized to 1 for g = 0. As there are D quartic melonic interactions (one for each cel?), one must introduce D 
intermediate matrix fields H c . The original tensor T and T can be integrated out to obtain: 


Z sym (g) 



S(H ) = \n d ~ 1 
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£Tr c [iTtf c ]-Tr c [lni?(tf)] , 


where Tr c is the trace over the index of color c, Trx> is the trace over all the indices, and: 
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where we denoted I c the identity matrix with indices of color c and l c = (g) c6C ;I c . 

In this paper we will discuss a quartic melonic model with only one interaction. In this case one uses only one 
intermediate matrix field and (dropping the index c) the partition function simplifies to: 
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Observe that in this form a tensor models in rank D looks very similar to a matrix model for alVxJV matrix. 
The only difference comes from the modified scaling factor TV D ~ 1 instead of TV for the action. This apparently 
innocuous modification leads to numerous consequences which we will be exploring in detail in this paper. 


3 Vacua of the theory 


The equations of motion of the model defined in Eq. © are: 
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1-gH 
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Due to the unitary invariance, Eq. © fixes only the eigenvalues A.; of H: upon diagonalization, the equations of 
motion become: 

A* - g = 0 for i= 1, - - -, TV, (8) 

l-g\i 


whose two solutions are: 
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We have a critical value of the coupling at g c = 1/2, where we meet a square root singularity and the two solutions 
merge. The solution a_ has a regular Taylor expansion around g = 0 (starting at linear order, a_ = g + 0(g 3 )), while 
the solution a + has a simple pole at the origin. 

Up to conjugation by an arbitrary unitary matrix, any solution of the equations of motion (i.e. a vacuum of the 
theory) is of the form: 


w->=(“o* „j_)=“*(o sms <io) 

where 1+ and I_ are the identity matrices of size TV+ x TV+ and respectively TV_ x TV_, with TV+ + TV_ = TV, otherwise 
arbitrary. There are TV + 1 possible solutions of this type. 

Observe that from the onset we have neglected the contribution of the measure [dH] to the equations of motion. 
In fact this contribution can profoundly alter the equations of motion and it is possible that the “vacua” we identified 
by just considering the action are nowhere near the true vacua of the theory. As we will see below, this is in fact 
what happens in the case of matrices D = 2. However, for tensors (D > 3) the contribution of the measure is 
subleading in 1/TV and the solutions in Eq. (fTOl) are true vacuum states in the large TV limit. 

Each TV(Ar +j jv_) is invariant under the residual symmetry group U(N + ) xU(N_): 

H {n+ ,n.) = H ( n+, N .) jJ t ) - 6 W) - v[/ - s W-) > (11) 
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and, being a representative in a conjugacy class of solutions of the classical equations of motions, is a saddle point 
of the action S(H). The partition function is as a sum over the saddle points H( N+>N _y. 
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In the rest of this paper we will discuss the path integral near a saddle point H (from here on we drop the 
indices ( N+. AT) in order to simplify the notation). 

Evaluating the action at a generic solution H , we obtain: 


S(H) = ~Tr[H 2 ] -Tr 


In (^)]- JV *(f- 1 "7) + Ar -(T- 1 "7) 


(13) 


Near g = 0, regularity of the free energy demands N+ = 0. However, at g = g c = 1/2, the free energy develops a 
singularity, and beyond that, solutions with generic N+ are possible. For g > g Cl it can be easily checked that the 
real part of (Tbil) is independent of N+ and 1V_. Therefore all these vacua contribute to the same leading exponential 
order, and in the broken phase we do not have at this level any reason to choose a particular value of N+. 

Close to a vacuum H 7 we make the change of variables H = H + M, splitting the field into a background H and 
a fluctuation field M. Observe that even if H is not hermitian (which happens if a ± are not real), the perturbation 
M remains hermitian, and the translation by H is a translation of the diagonal entries of H in the complex plane. 
Using the equations of motion we obtain: 
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therefore the effective quadratic term (defining the inverse covariance) for the fluctuation field is: 

i Tr [ M 2 ] - i Tr[HMHM\. 


Introducing the following definitions: 

M ++ = Pj_ MP+ Af + “ = P 4 .MP_, M = P_MP_ 
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(M ++ M + ~\ 
\M~ + M—J ’ 


and taking into account that P+ + P- = I, the quadratic part for the fluctuation field rewrites as: 

\{l - a 2 ) Tr[(Af ++ ) 2 ] + i(l - a 2 ) Tr[(M") 2 ] + (1 - a + a_) Tr[A^ + -M- + ] . 
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The inverse covariance is an operator in the vectors space of hermitian matrices with inner product ( M\\M 2 ) = 
TrfAfiA^]. Since: 

(P + ® P + ) | M) = \P+MP+) = | M ++ ) , (17) 

the inverse covariance is: 


(1 - a 2 )P+ ® P + + (1 - a?)P_ <8> P_ + (1 - a+a_)(P+ ® P_ + P_ ® P + ) , (18) 

and P+ ® P+, P_ ® P_ and P+ ® P_ + P_ ® P+ are (mutually orthogonal) projectors. The mass eigenvalues for the 
fluctuation field M are then read off as A+ = 1 - a 2 with degeneracy IV 2 , A_ = 1 - a 2 with degeneracy IV 2 , and 
1 - a+a_ = 0 with degeneracy N+N_. 

Defining e = 1 - 4t; 2 we have A+ = 2^/e/(-y/e t 1), and we distinguish the following situations: 
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• For small e > 0, the leading behavior of the mass eigenvalues of the fluctuation field is A ± ~ T2fye, and the only 
positive eigenvalue is A_. It follows that for g < g c only the solution AL = N,N+ = 0 is stable, and therefore 
we are in the U(N) symmetric phase, in agreement with the small g analysis. 

• At e = 0 (i.e. at g = g c ), all the mass eigenvalues are zero, signaling criticality. 

• For small e < 0, the leading behavior is A ± ~ 2|e| t i2y/jef, and both A+ and A_ have a positive real part. 
It follows that for negative e both vacua are stable and therefore we are in the broken phase, with residual 
symmetry U(N +) x£fyAL). In this case the off-diagonal fluctuations M+_ and M_+, lacking a quadratic term, 
have the interpretation of massless Goldstone modes. 


Around a vacuum state the effective theory for the fluctuation field writes: 
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The theory thus obtained looks quite complicated. We can considerably simplify it by getting rid of the off- 
diagonal blocks M + ~ and M~ + (which only appear at cubic order in the above action), exploiting the unitary 
invariance. The procedure we detail below comes to explicitly integrating out the Goldstone modes, and obtain the 
effective theory for the “radial” degrees of freedom. 


4 Integrating out the massless modes. 


Knowing that at the saddle point the eigenvalues split in two sets of N + and AL identical eigenvalues, it is convenient 
to consider the subgroup H - U{N+) x ZfyAL) c U{N) defined by: 
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for two arbitrary Hermitian matrices W+ and W- of sizes N+ x N+ and respectively N- x AL. 
We first show that any hermitian matrix H can be brought into block diagonal form: 

h - v [ H o l) V '’ 


( 20 ) 


( 21 ) 


for some V € U(N)/7i, i.e. for some V in the left coset of H in U(N). We subsequently perform the change of 
variables H -+ (H+,H~,V), and compute the Jacobian [ dH] = defy J)[dH+][dH-][dV]. 


Block diagonalization. Let us write H = UEU\ where E is diagonal and U e U{N) is a unitary transformation 
that diagonalizes H. Let us denote V € ZJ(7V)/H the representative in U(N)/H of the equivalence class of U, 
[U] = {UW\ W e %}. Then there exists a W eTL such thatQ: 

( -iW+ p iW + n \ 

0 e- iW E^ w ~y' (22) 

where E+ denotes the first N+ eigenvalues of H and E- denotes the subsequent AL eigenvalues. Defining H+ - 
e lW+ E+e~ lW+ and iL_ = e lW ~E-e ~ lW -, shows that H can be brought into block diagonal form by conjugation by 
some V eU(N)/T~L. 

1 For the sake of precision, U e U'(N) = U(N)/U( l)®fy because a multiplication of U to the right by a diagonal matrix 

diag(e l,il ,..., e z ^ N ) leaves H unaffected. Similarly, we should also quotient H ->■ H' = 'H/W(1)® JV . However, - U'(N)I'H'. 
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Jacobian. In order to compute the Jacobian J of the change of variables H -*■ (U + ,ff_, V), we will first show 
that the Jacobian does not depend on V and then compute it for V close to the identity. 

Denoting W+ and W- two Hermitian matrices, and B a complex matrix (of sizes N+ x N+, N- x N- and N+ x iV_, 
respectively), any U eU(N) close to the identity writes as: 


( W + B\ l 0 B\ ( W + 0 \ (W+ 0 \ 

U = eW W -)=e\ Bi °/el° W >+0{BW), el° W ~Kh , (23) 

hence the class of U in U(N)/'H can be represented by: 

?)• ( 24 > 

The infinitesimal variation of V is then: 

dv=i [d% ?)■ < 25 > 

and for any fixed W eU(N) close to the identity, W = I + iT (with T a hermitian matrix), we have: 

WV»I + iT + i^ +0(T 2 ,B\QB)^[d(WV)]*i[ d ° B i d *\ « [dV] , (26) 


hence the Jacobian does not depend on V. In fact, the invariance of the measure induced on the quotient space by 
the Haar measure on the group can be proved under quite general conditions (which in particular apply to our case 
as U(N) is a compact group), see for example [TT1 . 

We now follow the standard strategy. First, we write an infinitesimal variation of the matrix H as: 
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and use VV t = I => dV V' +VdV' = 0 to eliminate dV^. Then we use the invariance of the measure and evaluate 
the Jacobian at V = I: 
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Finally, writing dHij = Jij^kidXki, for X = (H+, H_, B, B^) and H = ( H ++ ,H ,H + + ) we read off the Jacobian: 


/I+ ®1+ 


J = 


dH 

dX 


0 


0 0 
0 I_ ® I_ 0 0 

0 0 i(I+®H--H + ®I-) 0 

0 0 0 -i(H_ 0 1+ -1_ 0 H + ) t 


(29) 


with determinant det( J) = det(iJ+ ® I_ - 1+ ® iJ_) 2 . We have thus shown that: 
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The integral over V is now completely decoupled and can be absorbed into a proper normalization. 

This formula generalizes the usual reduction to eigenvalues, which can be recovered by taking iV_ = 1, and then 
iteratively repeating the decomposition on the remaining N + x N+ block. This will iteratively generate the usual 
Vandermonde determinant A(A ) 2 = n,;<j(V;-Aj) 2 . Keeping instead N+ and V_ generic, but reducing to eigenvalues 
our formula (1301) . we write: 
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and of course the two sets of eigenvalues can be recombined as A = (A + , A ) to give the usual expression: 
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The effect of the scaling with N D 1 in Eq. © comes now into play. From Eq. (13211 we conclude that we 
can evaluate the partition function by a saddle point method without having to care about the Vandermonde 
determinant if D > 2. In fact log A (A) is of order TV 2 , whereas N D ~ l Y.iS(Ai) is of order N D . This is de facto what 
we have done when writing the saddle point equation © for the action alone (i.e. without measure contributions). 
Eq. © yields just two solutions for the eigenvalues at the saddle point. We can then directly use eq. (ED and 
perform the saddle point approximation, taking in each block one of the solutions i.e. A + = a+ and A~ = a_. In 
order to study the fluctuations around the saddle point one would then write Af = a± + fif and study the corrections 
coming from the diagonal fluctuations qif. 


5 The effective theory 

If one computes the partition function using the eigenvalue decomposition, the connection between the Feynman 
graphs of the theory and random triangulations is lost. In order to keep this connection explicit and find an 
interpretation of the effective theory around a non trivial vacuum in terms of random triangulations, one must start 
from Eq. m , and decompose the matrices H± into a background and a fluctuation H± = a± + M±. We thus obtain: 

Z(N+,N-) ~ e L ' ' ' 

x J~[dM+][dMS] |^det ((a+- a_)I + ® I_ + M + ® I_-1+® AL^J e~ N s , 

S eS =\ (1 - a 2 ) Tr[M 2 ] + \(1 - a?) Tr[M 2 ] - £ ^ Tr [M p + ] - £ ^ TV [Ml] , (33) 

Z 1 p> 3 P p> 3 P 

achieving our goal of eliminating the off-diagonal blocks in Eq. ED, at the price of introducing a non trivial 
determinant term. Note however that M ± in Eq. ( 1351 ) are not the same matrices as the M ++ and M in Eq. m, 
in particular they depend non-trivially on the off-diagonal blocks M + and M + in Eq. (1151) just like the eigenvalues 
of a matrix depend on its off-diagonal elements. 


5.1 Feynman graphs 

The integral in Eq. (1551) for the fluctuations around a saddle point: 


(a+ - a-) 2N+N ~ f [ dM + ][dM _] e ^[in(W- + ^(M + ® L - I+ ®M_))] 

-JV D -h(l-aJ)'&[Mj] + i(l- 0 ’) Tr[Mf]-E p > 3 ^ Ti[A^]-E p > 3 ^ Tr[Mf]l 

x e l j 5 

can be written more explicitly by expanding the logarithm: 
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(34) 


Observe that the terms coming from the non trivial determinant generate new quadratic terms. These terms 
could be included in the measure, but we prefer to treat them as bivalent vertices. 

The logarithm of the partition function is a sum over connected graphs. The graphs have three kinds of vertices: 
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• vertices coming from E p+ >3 Tr [Mf t ]. We denote v p+ the number of such vertices of coordination p+ and 
V+ the total number of such vertices, V+ = E P+ >3^ P +- These vertices have coordination at least 3. These 
vertices represent polygons in the + sector, which will glue together to form + surfaces. 

V— 

• vertices coming from Z p _>3 Tr We denote v p _ the number of such vertices of coordination and 

Vi the total number of such vertices, VI = Z P _> 3 ^ P _- These vertices have coordination at least 3. These 
vertices represent polygons in the - sector, which will glue together to form - surfaces. 

• vertices coming from q 1 ( a+ ~a-) q + +q - (% + + 9 ) T r [-^+ + ] Tr We denote v q+q _ the number of 

such vertices of coordination q+ and q- and V+- the total number of such vertices V+- = H q q+ + q I>Q v q+q _. These 
vertices have a joint coordination q+ + at least one, but can have coordination 0 in one of the sectors. These 
vertices represent a pair of a + and a - polygons connected at a nodal point. The + polygons (resp. the - 
polygons) will glue into + (respectively -) surfaces, and the + and - surfaces will be connected at a mixed 
vertex by a nodal point. 

We can represent the + and - vertices as solid vertices, and the +- vertices as two solid vertices (one for the + 
sector and one for the - sector) connected by a dashed edge, see Figure [l] Each graph represents a nodal surface 
(for definitions and properties of nodal surfaces see for example [42] ). The E+ edges connecting two M+ half edges 



Figure f: An example of a graph generated by the effective action in Eq. (l34l) . Here the two dashed edges represent 
two +- vertices of the type Tr [M 2 ] Tr [M 2 ] and Tr [M 2 ] Tr [M_], 


(which we call the + edges) bring a factor The edges connecting two M _ half edges (which we call 

the - edges) bring a factor • 

The graph is connected but, by erasing the dashed edges, it might disconnect into several connected components. 
Let us call C+ the number of the connected components of the graph made exclusively of the edges +, g+ the sum 
of the genera of these components and F+ the number of faces of these components (and similarly for the - 
components). Up to irrelevant combinatorial factors the amplitude of a graph is: 


( 1 \ ( 1 ^ N {D-1)V + Z p+ * 3 P+Vp + n (D- 1)V- Xp-x 

\ W°~i(l - a 2 ) / \ A r£>_1 (l - a 2 ) / 
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-Nf + N F ~ 


(35) 


Note that in the symmetric phase (i.e. g < g c and N + = 0) the amplitudes are positively definite, while in the broken 
phase (i.e. g > g c and N + > 0) they are complex. In the following we will not keep track of global minus signs for 
individual amplitudes. 

The various numbers of vertices, edges and faces are related by the following topological relations: 


q++q-> l 

2E + = Y, P+ v p + + Z 9+ u 9+9- > 
p+>3 q+,q->0 


q++q~>i 

2 E_ = Y P- v p- + Z 9-^9+ 9- > 

p~> 3 q+,q-> 0 


(36) 


U + + U + _ - E+ + F+ = 2 C + - 2 g + , 


V- + V+- -E.+F. = 2C- - 2g_ , 


(37) 


E + - V + - V+- + C + > 0 , 


E.-V.- U + _ + C. > 0 








E++ E--V+ -V- - V+- + 1 > 0 , 
V + . - C + - CL + 1 > 0 . 


(38) 


Equations (l36l) and (l37l) are usual graph relations applied to the + and - subgraphs: the first simply states that 
summing the numbers of half edges at every vertex amounts to counting twice the number of edges, while the second 
is simply the Euler characteristic formula. Equations (l38l) state that the first Betti number of a graph is always 
non-negative, in particular it is zero for trees and positive for graphs with loops (in fact it counts the number of loop 
edges, or excess edges). A generic connected graph, consisting of + and - components and dashed edges, is dual to 
a nodal surface with V+- nodal points. The total genus of a nodal surface is g+ +g- + ( V+--C+-C- + 1), i.e. the sum 
of the genera of the + components, the - components and the first Betti number of the abstract graph obtained by 
collapsing the + and - components to vertices connected by the dashed edges. Of course g++g-+(V+--C+-C-+l ) > 0, 
and the equality holds only for the case in which all the + and - component are planar, and the abstract graph is 
a tree. 

Defining the filling fractions n+ = ^ < 1 and n_ = ^ < 1, the amplitude of a graph writes: 

Sp + >3 P+Vp + E p _> 3 P- v P- 

_^^_ n F + n F- 

(1 - al) E +( 1 - a 2 ) E -(a + - ! o<9+ + «-K**- 

x n (D-1)V + -(D-1)E + +F +n (D-1)V + -(D-1)E + +F. 


which, using: 


(D - 1)V+ - (D - 1 )E + + F + = (D - 1)V + - (D - 1 )E + + E+-V+- V + _ + 2C + - 2 g + 
= -(D - 2)(E + - V + ) - E + _ + 2C+ - 2g+ , 


becomes: 


£p+>3 P+ V P+ Sp_>3 P— V p— 

a + a - F + F- 

- + ._>! - n + *n_ 

(1 - al) E +( 1 - a 2 _) E -{a + - 

x jy-(D-2)(E+-V+)-(D-2)(E--V-)+2-2(V + --C+-C- + l)-2g + -2g~ 


6 The phases of the model and their geometry 


(39) 


The main difference between the D = 2 case of matrices and the D > 3 case of tensors comes from the following 
theorem. 

Theorem 1. The number of graphs at fixed order in 1/N is: 

• finite for D > 3, 

• infinite for D = 2. 

Proof. Let us consider the scaling with TV of a term: 


-(D - 2 )(E + - V + ) - (D - 2)(E_ -V.) + 2- 2(T + _ -C+-C. + 1) - 2 g + - 2g_ . 


As the graph is connected V+- - C+ - C_ + 1 > 0, hence this scaling is bounded from above by: 


2-{D-2){E + -V + )-{D-2){E.-V.) . 


The V+ and V- vertices are at least three valent, hence 2E+ > 3b+ and 2E- > 3TA, and the scaling with TV of a term 
is bounded from above by: 

2 -£zl( E+ + E.). 

Finally, the number of dashed edges equals V+- and is itself bounded by E+ + E- + 1. Consequently, for D > 3, at each 
fixed scaling in 1/TV only connected graphs with at most a finite number of (+, -, and dashed) edges contribute. 
As the number of such graphs is finite the first assertion of the theorem follows. 
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For D - 2 the above bound does not constrain the number of graphs (it just states that the amplitude of any 
graph is bounded by TV 2 ). The scaling with TV of a graph becomes for D = 2: 


2 - 2 (V + _ - C+ -C. + 1 + g + + g.) , 

i.e. it is proportional with the total genus of the nodal surface dual to the graph. As the number of triangulated 
surfaces of fixed total genus is infinite, the second assertion of the theorem follows. 

□ 

Having estimated a bound for the scaling in TV of the graphs generated by the effective action Eq. (El , we can 
now address the question whether our expansion around the “vacua” in Eq. m is justified. 

The connected one point function (A Tr[M_]) , is a sum over graphs having an extra - univalent vertex 

representing the external field M__. The scaling with TV of such a graph is: 


2 )(E + - V + ) - (D - 2 ){E. -V.) + 2- 2(V + _ - C + - CL + 1) - 2 g + - 2 , (40) 


where VL includes the extra univalent - vertex. As 2 E- > 1 + 3(VL - 1). It follows that the scaling with TV of a 
connected one point graph is bounded by -(D - 2 )(E+ -V+)-(D- 2)[£’_ - (VI - 1)] < --AA ( E+ + E- + 1) hence: 



connected 


jo( i), 

\0(TV-^) , 


D = 2 
D> 3 ’ 


(41) 


that is the connected one point function starts at order 1 for D = 2, but is strictly smaller in 1/TV for D > 3. 
Similarly, we can look at the free energy: 

T = -jjD lnZ ("M = (y - ^ J ) + n- (y - In j j + , (42) 

where T contains the graphs for which we proved that an upper bound for the scaling is N~’ B 3 2 ( E++E ^. That is, 
the correction from the effective theory to the free energy computed at the saddle point is of the same order as the 
latter for D = 2, but it is subleading in 1/TV for D > 2. 


6.1 The matrix case D - 2 

It is now apparent why the D = 2 case is very different from the D > 3. Indeed, for D - 2 the critical behavior of 
the amplitude (1391) of an individual graph has no relevance: the free energy at fixed order in 1/TV is a sum of an 
infinite family of graphs, and acquires a critical behavior due to the criticality of this infinite sum. In particular, 
it is well known that this sum becomes critical at the critical constant g' c = ^=, which is smaller than the critical 

value g c - 1, at which the amplitude of an individual graph becomes critical. 

The fact that one still obtains an infinity of graphs at any fixed order in 1/TV shows that the contribution 
coming from the measure (the Vandermonde determinant) can not be neglected in this case. The “vacuum state” 
of a diagonal matrix with entries a+ and a_ which we obtained by neglecting the contributions of the measure is far 
from a genuine vacuum state of the theory, because, as we showed above, Tr[M_]) t d is of the same order 
inl/TVas(A T r[iV]} 

connected* 

6.2 The tensor case D > 3. 

For D > 3 one can derive the critical behavior of the free energy just by looking at the critical behavior of individual 
graphs, as the free energy is just a finite sum over such contributions. The diagonal vacua we identified are close 
to the genuine vacua of the theory as the connected one point function is subleading in 1/TV. 

The symmetric phase 

We first consider the case e = 1 - Ag 2 > 0. The only positive mass eigenvalue is A_, hence the unique stable vacuum 
corresponds to filling fractions n+ = 0,n_ = 1. In this case we have no + components, and no nodal points. In this 
phase the graphs are made of a single connected - component, and their amplitude is: 

Ep_>3 P-Vp- 

^ _ N 2-{D-2){E--V-)-2g- 

(1 -a 2 ) E - 
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From the onset we observe that the non planar graphs are subleading. Indeed, if a graph is non planar > 0, 
its amplitude is strictly smaller than the amplitude of a planar graph having the same number of edges £L and the 
same number of loop edges E- — V- + 1 which is planar. 

In the critical regime the amplitude of a planar graphs is: 

_J_ N 2-(D-2)(E.-V -) 5 (43) 

and, as all the V- vertices are at least three valent, we have V- < |E_, therefore this amplitude scales at most like: 

7V 2 (yeIV^) E . (44) 

The suppression in N of this amplitude can be offset by the enhancement at criticality. In the double scaling regime 
N -*■ oo, aT 1 = sJIN 3 fixed, arbitrary planar three-valent graphs contribute at the same order and dominate the 
free energy. 

This is to be contrasted with the double scaling limit in the symmetric phase for the tensor model which includes 
all the interactions, Eq. ©, discussed in [24]. In that case, in the double scaling regime either only trees decorated 
with self loop edges (tadpole edges) on the leafs contribute for D < 6, or arbitrary three valent graphs contribute 
for D > 6. The family of arbitrary three valent graphs is non summable, hence for D > 6 the double scaling series 
is not summable, like it is for matrices. For D < 6, the family of trees decorated with self loop edges on the leafs is 
summable. However, it is very likely that it corresponds to a branched polymer phase. 

For tensor models with only one interaction the dominant family in double scaling is summable, but it consists 
in all the planar three valent graphs. It is most likely that this double scaling regime corresponds to a Brownian 
sphere geometric phase of the model. 

The free energy is constructed by summing all the planar three-valent graphs with (up to combinatorial factors) 
weight a E ~ , therefore we expect to find a finite critical value a c at which graphs with an infinite number of edges (or 
vertices, i.e. triangles in the dual surface) dominate. This is the standard situation that one exploits to construct 
a continuum limit, but we will not push this further here. 


The broken phase 

When e = 1 - 4 g 2 < 0 both mass eigenvalues A+, A_ have a positive real part, hence we can have n+ > 0. 
regime the amplitude of a graph is proportional to: 

F+ F- 

_ U + n - _ jy-(D-2)(E + -V + )-( L D-2)(E--V-)+2-2(V+--C+-C-+l)-2g + -2g- 

^-E ++ E. + Y.lX + ,tl o' (9 + +<?-) «- 

As in the symmetric phase one can from the onset consider only the graphs such that g+ and g- are zero. 
Using the topological relations, the amplitude of a graph with planar + and - components is: 


In this 

(45) 


F F 

n + + n -~ /-^p+>3P+ v p+ T. p -23P+ v p- N -(D-2){E + -V + )-(D-2){E„-V-)+2-2(V + --C+-C- + l) 

r -3E++3E.^ /e JV 

\/e 


(46) 


Let us consider a graph Q and suppose that it possesses an internal V+ vertex of coordination p+ > 3. We compare 
the amplitude of this graph with the one of the graph Q' which is identical with Q, except that the p+ valent vertex 
is replaced by a 3 valent vertex and ap+-l valent vertex connected by a + edge. As Q' has one more + edge and 
one more + vertex with respect to Q, the scaling in N of the two amplitudes is the same. We have: 

A(g') = s/t +p+ ~ l ~ p+ A{g) => A(g') > A(g) . (47) 

\F 


It follows that the most singular graphs must have only three valent + and - vertices and their amplitude is: 

F F 

n + + n_~ /-3V++3V- pr-(D-2)(E + -V + )-(D-2)(E--V-)+2-2(V + --C+-C-+l) 

r-3E++3E- iV 


2+E+/3 2+E-/3 

= n. n_ ' 


N‘ 


/ /— D ~ 2 

(yiiv— 


.-(J5 + +B„) 


N' 


2(V + --C+-C-+1) 


(48) 


In the double scaling regime N -*■ oo, 


a 1 = s/eN°3 2 fixed, graphs with: 
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• only three valent V+ and V- vertices, 

• planar + and - components, 

• dashed edges that form a tree connecting the planar + and - components, 

contribute to the same order and dominate the free energy. Such graphs represents planar nodal surfaces, such that 
all the vertices which are not nodal points are three valent. 

Leading graphs in the double scaling regime in the symmetric and broken phases are represented in Figure [2j 







Figure 2: Examples of graphs dominating the double scaling regime. On the left a graph dominating the symmetric 
phase and on the right a graph dominating the broken phase. 


In the double scaling regime, the free energy in the broken phase is constructed by summing all the graphs 
detailed above, with a weight proportional to (n^ 3 a) E+ (n]/ 3 a) E ~ . Therefore it is effectively a double series in the 
parameters z+ = n 1 J 3 a and Z- = n\^ 3 a, which we expect to have a non-trivial domain of convergence in C 2 . For 
n+ = n_ = 1/2 and hence z+ = = z, this series becomes again a simple series in one variable, z. In this case we 

expect to find a finite radius of convergence. Note however that the coefficients of the power series are now complex, 
therefore one might have to tune z to a complex value in order to attain a continuum limit with a standard (in the 
sense of DT) geometrical interpretation. 
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